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ABSTRACT: We include finite top- and bottom-mass effects in the next-to-next-to-leading 
order parton shower (NNLOPS) event generator for inclusive Higgs boson production in 
gluon fusion based upon the POWHEG+MiNLO approach. Since fixed-order results for quark- 
mass effects only reach NLO accuracy, we add them to the NNLOPS generator at that 
accuracy. We explore uncertainties related to the unknown all-order logarithmic structure 
of bottom-mass effects by comparing the assumption of full exponentiation to no exponen¬ 
tiation at all. Phenomenological results showing the effects of finite quark-masses in the 
NNLOPS simulation are presented. These suggest that the aforementioned uncertainty is 
well contained within the envelope of plain renormalization and factorization scale uncer¬ 
tainties. 
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1 Introduction 

Since the 2012 LHC discovery of a new spin-zero particle [1, 2], all subsequent experimental 
analysis (see e.g. [3-11]) has revealed that its properties are very much consistent with 
those predicted for the Standard Model (SM) Higgs boson. The discovery of what appears 
to be the SM Higgs boson, without any new physics below the TeV scale to stabilize its 
mass, poses profound and challenging questions for the theoretical community. By the same 
token, LHC studies to quantify the level of agreement between SM Higgs boson predictions 
and experimental measurements, by increasing the precision in both, has further intensified. 

Thus far, tests and measurements of Higgs boson properties, in almost all cases, have 
focused on rather inclusive observables, such as production cross sections and branching 
ratios. However, even for fully inclusive quantities experimental analysis do often classify 
events in categories depending on the number of accompanying jets, and, at times, jet-vetos 
are used to enhance the signal-to-background ratio. In general these measurements require 
the use of fully exclusive Monte Carlo generators, possibly in combination with dedicated, 
theoretical calculations. 
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For gluon-fusion Higgs production, first differential measurements have been published 
in the second half of 2014 by the ATLAS collaboration [12, 13]. These measurements are at 
present statistically limited, however, forthcoming LHC Run-II data will lead to substantial 
improvements, and by the end of Run-II the measurement uncertainties are estimated to be 
substantially reduced. Full exploitation of these data, in particular the degree to which we 
are able to resolve or exclude new physics in the Higgs sector, requires that all theoretical 
tools, including Monte Carlo simulations, be as precise as possible. 

Next-to-leading order calculations matched to parton showers (NLOPS) [14-16] are 
now the standard for fully exclusive Monte Carlo predictions and are indispensable tools 
for the LHC experimental collaborations. 

Recently, NNLOPS generators for Higgs [17, 18] and Drell Yan [19, 20] production have 
also appeared. In particular, the generator of ref.[17], based upon the MiNLO method [21], is 
particularly attractive since it achieves NNLO accuracy without recourse to an unphysical 
partitioning of phase space. This generator is based upon the large mt approximation of 
the Higgs coupling to gluons. On the other hand, it is known that finite quark-mass effects 
are important especially in the case of production in association to energetic jets, since the 
radiation resolves the internal structure of the gluon-Higgs effective coupling. 

Finite quark-mass effects to the total cross section are known exactly at NLO (O (a|)) 
in QCD perturbation theory [22-27]. Finite quark-mass effects on high transverse momen¬ 
tum Higgs production were first calculated at LO (O (ct|)) in refs. [28, 29] and, as an 
expansion in 1 /mt, at NLO for finite top-quark-mass effects by Harlander et al. [30]. Some 
of the former finite quark-mass corrections have been implemented in public fixed order 
computer codes e.g. MCFM [31] FehiPro [32, 33] and SusHi [34], 

Finite quark-mass effects have been included in NLOPS Monte Carlo simulations of 
Higgs boson production in POWHEG [35] and MC0NL0 version 4.08 onwards [14], Recently 
the Sherpa collaboration have introduced a treatment of finite top-mass effects [36] in the 
context of their multi-jet NLOPS merging scheme MEPS@NLO. Differences between the 
predictions of the POWHEG program of ref. [35] and MC0NL0, in particular regarding their 
response to the inclusion of finite 6-quark-mass effects, have been a source of discussion and 
have stimulated further analytic work on the treatment of bottom-mass effects in resummed 
calculations. 

The earliest analytic resummation work to include mass effects was that of Mantler 
and Wiesemann at LO+NLL accuracy [37], work which was later extended to include also 
MSSM effects [38]. An NNLL+NNLO resummation of the transverse momentum spectrum 
of the boson in the larg e-mt limit was subsequently augmented with finite t- and 6-quark- 
mass effects at the NLL+NLO level by Grazzini and Sargsyan [39]. Shortly following this 
work Banfi et al. presented an extension of their NNLL+NNLO large-m^ computation of 
the efficiency for Higgs boson production in the presence of a jet veto in ref. [40]. All 

of these analytic resummation computations took a different approach to the inclusion of 
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finite 6-quark-mass effects, in particular, to their handling of the enhanced ~ 4j- log 2 ^ 
terms in the region m f, < px < rn ff. 

In this work we address the extension of the NNLOPS event generator for Higgs boson 
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production in ref. [17] to include finite top- and bottom-quark-mass effects. In section two 
we present our theoretical rationale and method for doing this, with each subsection describ¬ 
ing the implementation of finite quark-mass effects in a different layer of the event generator. 
Sec. 2.1 concerns the required modifications at the level of the underlying HJ NLOPS gen¬ 
erator [41]. Sec. 2.2 enters into the discussion on the role of finite quark-mass effects and 
resummation, describing their implementation (or not) in the MiNLO Sudakov form factor 
and the theoretical justification for our approach. Sec. 2.3 describes how the finite quark- 
mass effects are accounted for in the NNLO predictions required for the reweighting stage. 
In Sec. 3 we present a selection of numerical results obtained with the new NNLOPS gen¬ 
erator using various options to explore theoretical uncertainties. Conclusions are drawn in 
Sec. 4. Some technical details relevant to our new implementations are illustrated in the 
Appendix. 

2 Method 

In ref. [17] a generator that is NNLO accurate and includes parton shower effects (NNLOPS 
accurate from now on) was presented for inclusive Higgs production. This generator is 
based on the large nit effective theory, in which the Higgs boson emitted from a top-loop is 
approximated by an effective tree-level coupling of the Higgs to gluons and the interaction 
of the Higgs boson to lighter quarks (including the bottom quark) is neglected. 

Our goal here is to correct this NNLOPS Higgs generator in such a way that at least 
the most important mass effects are included. The NNLOPS generator relies upon three 
main components: 

• the POWHEG HJ generator of ref. [41]; 

• the MiNLO procedure on HJ discussed in ref. [21]; 

• the NNLO fixed order calculation of ref. [39]; 

• the reweighting procedure described in ref. [17]. 

In the following we describe how, and to what extent, we include mass effects in the first 
three items of this list. The reweighting procedure remains obviously the same. 

2.1 Mass effects in the HJ generator 

An NLO calculation (i.e. of order a^) of Higgs plus jet production including the finite 
quark-mass effects is not available at the moment, and we will thus resort to the following 
approximations. We begin by considering the quark-mass corrections due to the top loop 
only. We multiply all components (Born, virtual and real) of the infinite top mass approx¬ 
imation formulae by the ratio of the matrix elements for the production of a Higgs boson 
in association with a light parton, including exact one loop top mass dependence, divided 
by the same matrix element in the infinite top mass approximation. These matrix elements 
are all evaluated at the Born kinematics for the Born and virtual contribution, and at the 
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underlying Born kinematics for the real one. 1 This procedure guarantees that at the Born 
level the cross section has the full top mass dependence. The real emission cross section in 
the soft and collinear limits, has also fully corrected top mass dependence, and the same 
holds for the part of the virtual corrections that arises from virtual gluons of energy much 
below the top mass. Thus, the contributions that we miss have to do either with real emis¬ 
sions with widely separated jets, or virtual corrections where the gluon energy is not small 
with respect to the top mass. These corrections are in fact of two-loop level, and are not 
yet available. 

The inclusion of bottom mass effects is more delicate. In this case, in fact, there is 
still factorization for light parton emissions at scale below the bottom mass. However, this 
scale is now much smaller than the Higgs mass, and such factorization may turn out not 
to be useful in most of the range of light parton emission, or of light parton virtualities in 
one-loop corrections. We thus consider two options in this case. In the first one we ignore 
this problem, and include the effect of the bottom quark in our rescaling factor. In the 
second option, we apply the top mass correction to all components of the cross section, but 
include the bottom mass effect only in the Born case. 

The matrix elements for H +jet production at order a]? and with full quark-mass de¬ 
pendence have been computed in ref. [28, 29], and are implemented in the POWHEG BOX in 
the generator gg_H_quark-mass-eff ects [35]. We have used the code of ref. [35] for our 
purposes. 

We perform our reweighting at the level of the calculation of the matrix elements. The 
following reweighting factors are computed together with the Born term: 


where 


K t b = 


Btb(®B,m t ,m b ) 


B 


K, = 


inf 


B t ($B,r n t ) 

-®inf 


B[nf — ) Wit i ™b) |t?it=oo,m(,=0 • 


( 2 . 1 ) 

( 2 . 2 ) 


These reweighting factors are stored, and the Born matrix elements are multiplied by K t j,. 
The real and virtual corrections are multiplied by in our first option, and by Kt in our 
second option. 2 

We notice that with this procedure, the Born matrix element is computed with the exact 
quark-mass dependence. At variance with other methods, however, also the higher order 
terms are reweighted point by point in the phase space, either using the same reweighting 
factor, or a factor including only top mass effects. Our motivation to do so is that this is 
correct at least for the Born term and for the terms proportional to it. 


2.2 Mass effects in the MiNLQ procedure. 

MiNLO [42] is in essence and extension of the well-known CKKW procedure [43] at the 
NLO level. One first associates a most likely branching history to the kinematic structure 

1 The underlying Born kinematics of a given real emission kinematics configuration is obtained by an 
n + 1 to n mapping procedure (specified in details in ref. [15, 16]) that is such that in the singular collinear 
limits corresponds to the Born kinematics of the factorized cross section. 

2 The Born matrix elements that are used to compute the real counterterms, the collinear remnants and 
the soft collinear terms are always rescaled in the same way as the real and virtual terms. 
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of the event. The hardness of the branchings are then used to set the factorization and 
renormalization scales. Furthermore, appropriate Sudakov form factors are supplied to the 
process, in full analogy with shower algorithms. Care has to be taken to subtract NLO terms 
arising from the expansion of the Sudakov form factor in order to maintain NLO accuracy 
of the cross section. The Sudakov form factors guarantee that the MiNLO improved cross 
section can be integrated down to vanishing momenta of the associated jets. Furthermore, 
in ref. [21], it was shown that, in the case of the production of a boson in association with 
a jet, by a suitable refinement of the procedure, one could achieve NLO (i.e. a ]?) accuracy 
for fully integrated quantities, i.e. even without requiring the presence of the jet. 

The MiNLO Sudakov form factor of ref. [21] is computed assuming that the main pro¬ 
duction vertex is pointlike, which is the case in the infinite top-mass approximation. This is 
a justified assumption if we only consider the top loop in the Higgs production vertex. Soft 
gluons are in fact characterized by energies below the Higgs mass. The flow of their momen¬ 
tum through the top loop is thus expected to have a very limited effect, in view of the large 
value of the top mass. On the other hand, also the bottom contributes to the Higgs coupling, 
with the interference of the bottom and top loop affecting the total cross section by about 
-7%. Furthermore, it was noticed that the POWHEG generator gg_H_quark-mass-eff ects 
yielded a transverse momentum distribution for the Higgs differing substantially with the 
analogous MCGNLO simulation, and with several analytic treatments of the Higgs transverse 
momentum spectrum. These differences are easily traced back to the fact that in POWHEG 
the Sudakov form factor is computed including the effects of the bottom finite mass. It is 
in fact given by the formula 


exp 


' P t >Pt 


B($ b) 


(W r 


rad 


(2.3) 


where R is the real emission cross section and B is the Born term. On the other hand, 
in standard resummation formulae the real cross section is replaced by its Altarelli-Parisi 
approximation [44], This assumes that the underlying Born cross section is not affected by 
virtualities of the incoming partons smaller than the Higgs mass. 

This issue was dealt with in different ways in the literature. In ref. [39] it was assumed 
that the resummation scale in the case of bottom mediated Higgs production should have 
been set to the mass of the bottom rather than to the mass of the Higgs, in order to 
satisfy the basic resummation assumption that the soft gluons should be softer than the 
quark-mass. On the other hand, in ref. [40], it was shown by a detailed analytical study, 
that at order the double logarithmic structure of the gluon emission in Higgs production 
with quark-mass effects remains the same as in the infinite mass limit. There are however 
differences at the single logarithmic level. The all order logarithmic structure, accounting 
for the complex three-scales problem induced by the bottom contribution, is at present not 
known even at the double logarithmic level. It is thus clear that this fact will introduce an 
uncertainty in our current prediction that must be modeled in some way. In the present 
work, we do the following. We consider two possible approaches to the computation of the 
Higgs Sudakov form factor. In the first approach, that we will label as MEMB (standing 
for Matrix Elements mj), we keep the same MiNLO Sudakov that we use in the traditional 
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HJ-MiNLO generator, thus assuming that no large logarithms arise at order higher than a ^ 
due to the inclusion of the bottom mass effects. In the second approach, that we call 
RMB (standing for Resummed mj), we assume that all what is present at order due to 
bottom mass effects exponentiates. In order to do this, for every Higgs production event, 
we multiply the MiNLO Sudakov form factor by 


Since 


A mb (pt,H,yH) = exp 




r Rt 

'Pt,H>Pt,H Bt 


dh>: 


rad 


B tb B t ) rad 

(2.4) 

R m t =oo d$ , 
d rad 5 

iDmt =oo 

(2.5) 


we notice that if we use the Sudakov form factor of eq. (2.3), evaluated in the infinite 
quark-mass limit, instead of the MiNLO Sudakov, the A mb factor turns it into the Sudakov 
of eq. (2.3) including bottom mass effects. Thus, this second procedure should yield results 
that are closer to those of the gg_H_quark-mass-ef f ects generator. 


2.3 Quark-mass effects in HNNLO 

Quark-mass effects have been implemented recently at NLO in the version 2 of the HNNLO [39] 
code. For the NNLO predictions used in this paper, we set approximV which corresponds 
to having exact top and bottom mass dependence at NLO and no approximate mass effects 
at NNLO. 


3 Numerical study 

We begin by considering the HJ-MiNLO generator, with the quark-mass effects included at 
different levels. We consider Higgs production at the 8 TeV LHC. We set rn H = 125.5 
GeV, and use in the following MSTW2008NNLO parton distribution functions [45]. For 
the heavy quark-masses we use nit = 172.5 GeV and m b = 3.38 GeV. The choice of the 
bottom mass deserves some discussion. 

3.1 The choice for the bottom mass. 

In our computation, quark-mass renormalization plays no role. Thus, we have freedom in 
our choice of the quark-mass scheme that we should use, either an on-shell scheme, or the 
MS or DR scheme. The one-loop relation among the masses in the different schemes is as 
follows: 

m m (Q) = JUpole x 

^Dr(Q) = "Vie X 

Choosing a s = 0.115 and m po i e = 4.75 GeV, for Q = 125.5 GeV we get 3.38 and 3.32 GeV 
for the MS and DR schemes respectively. 


1-— ( log -F 4/3 ] + 0(a 2 ) 

7r \ m 1 


1 - — log 


''pole 

Q 2 

*pole 


+ 5/3 + 0(a 2 ) 


(3.1) 

(3.2) 
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One could argue that, since the bottom appears in a loop involving momenta of the 
order of the Higgs mass, it should be more appropriate to use the bottom MS mass evaluated 
at the Higgs mass scale. The following observation also supports this view. We compute 
the Higgs total cross section 0 T O I NLO ( us i n g the gg_H_quark-mass-effects generator) 
including only the top contribution at LO and NLO level. We then compute the same cross 
sections, cr^ 0 ^ NL0 , including also the bottom loop, in the MS and on-shell scheme. The 
results for the ratios are shown in table 1. We see that using the MS value at the LO level 



on-shell scheme 

MS scheme 

/yLO /^.LO 
a tb / a t 

0.89 

0.93 

_NLO/_NLO 
a tb / 

0.93 

0.94 


Table 1. Effect of the mass schemes on the inclusive Higgs cross section at LO and NLO level. 

yields a mass effect that is closer to the one obtained at NLO. In the following, we will thus 
use the MS value. On the other hand, we have verified that changing the scheme for the 
top mass changes the cross section only at the per-mille level, both at LO and NLO. 

3.2 Mass effects in matrix elements 

For the following plots, we only consider results at the level of the Les Houches output, i.e. 
no shower effects are included. 

First of all, we would like to assess the difference between the two variants for the im¬ 
plementation of mass corrections in the matrix elements, discussed in Sec. 2.1, i.e. whether 
the bottom-mass correction should be only applied to the Born term, or to the full NLO 
cross-section. In figure 1 we compare the two procedures for both the Higgs rapidity dis- 



Uh Pt,h [GeV] 

Figure 1. Rapidity and transverse momentum distributions for Higgs production at the 8 TeV 
LHC. The two lines represent the result from the HJ-MiNLO generator with bottom and top mass 
effects fully included in the matrix elements (red, dashed line), and the same generator without 
bottom mass effects in the matrix elements for the NLO corrections (green, solid). 

tribution and the Higgs transverse momentum. As we can see, the differences are at most 
2%. We will thus, as our default option, include also the bottom quark in the rescaling 
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factor for the NLO corrections. The other option is still provided in the public code. It 
can be used to assign a systematic error associated to using an approximated treatment of 
bottom-mass effects at NLO in HJ. 

We now compare the plain HJ-MiNLO generator (without any quark-mass effects) to 
our new generator with mass effects included in the matrix elements only (and no mass 
effects included in the MiNLO Sudakov form factor), corresponding to our MEMB option. 
These are displayed in fig. 2 for the rapidity and transverse momentum distribution of the 




Figure 2. Rapidity and transverse momentum distributions for Higgs production at the 8 TeV 
LHC. The three lines represent the result from the plain HJ-MiNLO generator, and the result from 
the same generator improved with either the top mass effect alone, or with both top and bottom 
mass effects. 

Higgs boson. From the rapidity distribution, we notice that the inclusion of the top mass 
alone amounts to a constant increase of about 5%. On the other hand, the inclusion of the 
bottom-mass effects decreases the cross section by about the same amount, thus yielding 
a net negligible effect on the total cross section. A very similar pattern is observed in the 
NLO total Higgs production inclusive cross section. For example, ref. [39] quotes a 6 % 
increase due to the top mass, and a - 1 % effect when also bottom effects are included. 

The similarity between the HJ-MiNLO and the fixed-order inclusive NLO results for the 
total cross section is not trivial. In fact, the HJ matrix elements are those for Higgs pro¬ 
duction in association with one parton at NLO, and with two partons at LO, while the 
fixed-order inclusive NLO case involves the matrix elements for the production of the Higgs 
without accompanying partons at NLO, plus the production of the Higgs with one extra 
parton at LO. In general, the MiNLO procedure guarantees that even when the extra par- 
tons in the HJ generator are integrated out, the same accuracy of the fixed order, inclusive 
NLO calculation is achieved. In the present case, since the mass corrections in HJ-MiNLO 
are included only in an approximate way, mass effects in the inclusive cross section do not 
have full NLO accuracy. As we can see, however, the same pattern of the NLO calculation 
is recovered. From the transverse momentum spectrum we notice that, as expected, the 
difference between including only the top, or both the top and the bottom quark is concen¬ 
trated at low transverse momenta, below pt,h = 50 GeV. At zero transverse momentum 
this difference is largest and amounts to about 12 %. 
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3.3 Resummed mass effects 


We still consider results at the Les Houches level, and compare the plain HJ-MiNLO generator 
(without any quark-mass effects) to our new generator with mass effects included in the 
matrix elements, and also including the correction to the MiNLO Sudakov given in eq. (2.5), 
corresponding to our RMB option. 

Before doing so, we remind the reader that we expect these results to match to some 
extent those of the gg_H_quark-mass-eff ects generator. We thus begin by showing in 
fig. 3 results obtained with the gg_H_quark-mass-effects, when only the top loop is 




Pt,h [GeV] 


Figure 3. Rapidity and transverse momentum distributions for Higgs production at the 8 TeV 
LHC. The three lines represent the result from the gg_H_quark-mass-effects generator, when 
only the top loop is included, and when the bottom loop is added in the MS and in the on-shell 
scheme. 

considered, in comparison to the case when also the bottom loop is included. In this last 
case, we display results in both schemes, the onshell and MS scheme. As explained before, 
in our opinion the latter choice is preferable in this case with respect to the on-shell scheme. 
We see that, in the MS scheme, in the rapidity distribution the difference is very small. 
This is due to the fact that scheme compensation takes place in inclusive quantities, i.e. 
the explicit modification of the virtual contribution due to the scheme change compensates 
the variation induced by the change in the mass parameter. This compensation does not 
take place in the transverse momentum distribution, that can be considered a leading order 
quantity, mostly affected directly by the change in the bottom mass. Furthermore, there is 
an explicit m^ dependence in the Sudakov exponent, where no scheme compensation occurs. 
In the onshell scheme we find more pronounced differences, both at low and high transverse 
momenta. 

We now show in fig. 4 the HJ-MiNLO results in our RMB approach (red line) compared 
to including only top-mass effects (green line). As anticipated, the bottom-mass effect is 
now very similar to the one displayed in fig. 3. 

3.4 NNLOPS results 

We now turn to our full NNLOPS results. The NNLOPS events are showed with PYTHIA 
(v.6.4.25)) with the 2011 Perugia tune (PYTUNE(350)). Results include full hadronization 
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Pt,h [GeV] 


Figure 4. Rapidity and transverse momentum distributions for Higgs production at the 8 TeV 
LHC. The two lines represent the result from the HJ-MiNLO generator, when only the top loop is 
included, and when the bottom loop is added (in the MS scheme) in our RMB option. 


and underlying event effects. 

We start from the Higgs rapidity and 
displayed in fig. 5. 



transverse momentum distributions, that are 



Pt,h [GeV] 


Figure 5. Rapidity and transverse momentum distributions for Higgs production at the 8 TeV 
LHC. The two lines represent the result from the NNLOPS generator, in the large to* approximation 
(red), including only top-mass effects (green), and including top and bottom mass effects in the 
MEMB scheme (blue) or in the RMB scheme (magenta). 

We see that, for the rapidity distribution the MEMB and RMB results are very close to 
each other, and in fact very close to the large mt result. On the contrary, results including 
only top loops are about 6 % larger. This is in fact the result that one obtains at pure 
NNLO level. In the transverse momentum distribution on the other hand, we observe a 
difference, of up to 5%, between the RMB and MEMB schemes. 

In fig. 6 we show the four predictions for the leading jet integrated cross-section (left) 
and jet-veto efficiency (right). We notice that the difference between the MEMB and RMB 
schemes is very small (of the order of 1 - 2 %) for values of the transverse momentum of the 
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Figure 6. As in Fig. 5 but for the leading jet integrated cross-section (left) and jet-veto efficiency 
(right). 


order of 25-30 GeV, the region of interest in Higgs studies involving a jet-veto in ATLAS 
and CMS. A similar conclusion was found in ref. [46] where a resummation for the jet-veto 
was presenting including various matching procedure to estimate the uncertainty due to 
the top-bottom interference contributions. On the other hand this difference rises at lower 
value of the jet-veto cut, and reaches about 5% at PT,cut = 5 GeV. 

Finally, in fig. 7 we show the effect of parton shower and hadronization for both the 



Figure 7. Comparison of the result after shower, hadronization and inclusion of the underlying 
event (red), after shower only (green), and at the Les Houches level (blue) for the leading jet 
integrated cross section. On the left, the mt = oo result is shown, and on the right the RMB result 
is reported. 


mt = oo case and the mt,mb, RMB result. As we can see, the parton shower effect 
has a noticeable impact, while hadronization and the underlying event do not affect this 
observable sensibly. In this, we see no difference between the mt = oo and the mt,mb, 
RMB case. 
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4 Conclusions 


In this work, we have included finite quark-mass effects in our NNLOPS generator for Higgs 
production. As far as mass-corrections are concerned, the accuracy of our generator is 
that of currently available fixed-order calculations, i.e. NLO accuracy. In our procedure 
we have considered two possible sources of uncertainty. The first one, has to do with the 
fact that we rescale the full matrix elements with the mass effect that are computed at the 
level of Higgs production in association with one jet. This procedure may be considered 
safe as far as top mass-effects are concerned, where one may argue that dominant NLO 
corrections, involving momenta softer than the Higgs mass, will not affect the top-quark 
loop. This assumption is unjustified when the bottom loop is involved. We thus consider 
two alternatives: we either apply the mass-correction only to the Born contribution, or 
to the full NLO matrix element. We find that numerically, the two procedures lead to 
differences of the order of 2%. A second source of uncertainty has to do with the all¬ 
order inclusion of contributions enhanced by powers of log(mj/rojf). Little is known about 
the all-order structure of these logarithms. We thus consider two options: either we do 
not include them at all, or we fully exponentiate them in the MiNLO Sudakov form factor. 
The first approach is similar in spirit to applying a resummation procedure by assuming 
that soft gluon momenta do not affect the quark loops. The second approach is analogous 
to what is implemented in the gg_H_quark-mass-effects generator, where the POWHEG 
procedure leads to full exponentiation of the real emission cross-section including mass 
effects. In this case we find differences in the low transverse momentum region of the order 
of 5% starting from the low momentum region and changing direction at high transverse 
momenta. On the other hand, when considering jet-veto distributions for typical values of 
the transverse momentum cuts used by ATLAS and CMS (pT,cut = 25-30 GeV), we find 
only 1-2% differences. 

Unlike in ref. [17] we have not reported here a study of scale variation, since in view of 
the smallness of the quark-mass effects, we do not expect to get results sensibly different 
from those of reference [17]. 

The code for this generator has been made available in an update of the HJ package in 

the P0WHEG-B0X-V2. 
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A Technical details about 6-mass effects in the MiNLO Sudakov 

The inclusion of bottom-mass effects in the MiNLO Sudakov, following eq. (2.5), requires 
for each point an additional integration over a ratio a matrix elements, for fixed y#, over 
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the phase space region where the Higgs transverse momentum is larger than the Higgs 
transverse momentum of the event, px- Since such integration is performed on the fly for 
each point, it is crucial that the phase space is sampled efficiently. Furthermore, we will use 
an approximate expression of the matrix elements appearing in eq.(2.5). In Appendix A.l 
we give all details about the phase space evaluation, while in Appendix A.2 we illustrate 
the approximations done on the matrix elements in eq. (2.5). 

A.l Phase space 

Here we study the phase space integration of the matrix element for the production of a 
Higgs boson plus one light parton, under the constraint that the Higgs boson has fixed 
rapidity yu, and that its transverse momentum is larger than a given lower limit px- 

We denote by run the Higgs mass, kn = (E,k) its four momentum, by l g = (1,1) the 
light-parton momentum and define q = ku + lg • The three-vectors l and k are defined in 
the partonic center of mass (CM) frame. We have also defined k = k and 1 = 1. The two 
body phase space in the rest frame of the Higgs-light-parton system is then given by: 


d $2 = 9/(27r ) 3 27n H(g - lg) 2 ~ ™h) 



(A.l) 


where s = q 2 and where we have used the property k = —l. The ± sum refers to the two 
possible signs of the li integral (the suffix l stands here for longitudinal). 


The full phase space integral, to be multiplied by parton distribution functions and 
partonic cross section, including the constraint that the Higgs transverse momentum is 
larger than px, is then 



where y H is the Higgs rapidity in the laboratory frame and y^ m the Higgs rapidity in the 
CM frame. The above expression is equivalent to 


h k 



cm 


H 


I y H — ycm\)d(k 2 — kx)9(kx — Px) ■ (A.3) 


The Higgs rapidity in the CM frame can be written as 



(A.l) 
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where E. k are the energy and momentum of the Higgs in the partonic CM frame: 


_ s + mjr 1 + z 

E = — „ ^ = vtih- 


2\/s 


2 v/i’ 


k = 


s — m 


H 




= rn H 


1 - z 

Vi’ 


and s = m 2 H /z. Furthermore we define 


(A.5) 



s my 

5 = 5z ' 

(A.6) 

In order to solve for k'^ 

in the delta function we get 



\Jk 2 — = E tanh \y H - y cm \ , 

(A.7) 

and 

k t = k 2 - E 2 tanh 2 (y H - y cm ) . 

(A.8) 

Since we must have px < kx < k, from the theta functions we get the constraints 



p 2 T <k 2 - E 2 tanh 2 (y H - y cm ) < k 2 , 

(A-9) 

that is to say 

k 2 -Pt> E 2 tanh 2 {y H - y cm ), 

(A.10) 


or 


Ay = atanh 


k 2 - PT 2 _ 1 . (E+yJk 2 -^ 


E 2 


= 2 l0g 


> I?/# - Vcm | • 


The jacobian for the delta function integration is 

2(m 2 H + A;|.) \Jk 2 — kj. 


dy, 


H 

cm 


dkj* 


-l 


= -^(F ; 2 - k 2 + kl)^k 2 -k 2 T = 


E 


Our phase space integral becomes 

f , f-h l °sr 1 E 2 -k 2 + kl 

/ dr / dy cm --^7=- 

J J i logr 8vr 

which is equivalent to 


6 »(Ay - ly# - y C m|) 6 / (A : 2 - Pt) 


-1 


' TO 


, p logl / r , z m 2 +fc 2 

dr / dy cm -——- 

Ji logr 4vrm^(l + z) 


H "h AT #(F 2 - T 2 tanh 2 (i/H - y cm ) - p|), 


(A.ll) 


(A.12) 


(A.13) 


(A.14) 


with to = (jpt + Mt) 2 /S, with Mt — \jvr T m #- This limit arises from the requirement: 

k=s ~^ L>PT ' (A - 15) 

that implies 

s — 2 y/spr — m 2 H > 0 =>■ (y/s — {jp T + Mr)) (\/i — (pr — Mt)) > 0 , (A.16) 


- 14 - 



which in turn implies 


The theta function 


s > (pt + M t ) z 


0(k 2 - E 2 tanh 2 (y H - y cm ) - p\) 

also implies some restrictions on y cm , that we now work out. It can be written as 

I k 2 — p\ 


atanh 


E 2 


> |VH - Vc m|, 


or 


y H + A > y cm > y H - A, 


with 


A = atanh 


k 2 — Pt 


E 2 


Combined with the integration limits for y cm , this yields 

max log ^,y H - a) < y cm < min Q log ^ , y H + A 


The range is not empty if 


that can be combined into 


1 . 1 

-xlog- < y H + A, 

2 r 

VH ~ A < - log-, 
z r 


A > “log - + \y H \- 

2 r 


Taking the hyperbolic tangent of both sides we get 


k 2 — p^, — 2 7" — g 

.E 2 re 2 \vh\-\-\ t + e~ 2 \ yH \ 

Let us define r/ = we have the inequality 


A ; 2 -" 2 


E 2 


P T T ~y 
^ , ) 
T + p 


or 


( tS — m 2 H ) 2 — 4 Sp^T t — r) 


Y (tS + m 2 H ) 2 

that is easily solved to yield either r < rj, or 


t + rj 


t > n = 


(7/5 — m^) y/ rjSMj, + V^Pt 


5(t/5 - M 2 ) 


(A.17) 

(A.18) 

(A.19) 
(A.20) 

(A.21) 

(A.22) 


(A.23) 


(A.24) 


(A.25) 


(A.26) 


(A.27) 
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or 


— (pS — rn 2 H )yJr]SMj, + pSp, 


t < r 2 =- 7 — -7T - . (A.28) 

S(pS — M^) V ; 

Setting p = x 2 , and using S = E 2 and m 2 H = Mj, — p^, we get (using an algebraic 
manipulation code) 

x(EMtx + Pt + Mj.) 


T\ = 


T2 = - 


E(Ex + Mt ) 
x(EMtx — Pt + Mt) 


< 0 , 


(A.29) 


(A.30) 


E(Ex + Mt) 

so that only n is retained, and we need to satisfy the conditions: either r > To and r < 77 , 
or t > To and t > T\. We also find: 

n - „ = - J Ex ~ PT -Jt T } {E l\ PT ~ Mt) = -x (x - CT) £ A + / T ~ Mr . (a. 31 ) 


E(Ex — Mt) 


Ex — Mt 


and 


n - t 0 


M t { x -^) 2 
Ex — Mt 


(A.32) 


We notice that we must have Ex — Mt > 0 in all cases, otherwise the Higgs with transverse 
mass Mt and rapidity yn would be inconsistent with the total incoming energy. Hence 
t1 > tq, and we get only two cases: 


1. To < p: in this case the right hand side of eq. (A.31) is negative. In fact, E^/tq — Mt is 
obviously positive, and it remain positive if we replace ^/tq with x, since x > ^/tq. So 
the last factor on the right hand side of eq. (A.31) is positive, and the right hand side 
is negative. Therefore ti < r/, and the range of integration in t is To < t < 1, because 
t > rj also implies t > t\. 


2. tq > p: in this case eq. (A.32) guarantees that ti > to, and we have T\ > tq > p. The 
range of integration is ti < t < 1. The condition T\ < 1 is the remaining requirement of 
consistency for the assigned yn and pr values. 


A.2 Approximate matrix elements in the MiNLO Sudakov 

Here we record the approximate matrix elements used for computing our finite quark-mass 
corrections to the MiNLO Sudakov form factor. For what concerns quark-mass effects in 
the partonic Higgs boson-plus-one parton sub-processes, all other components entering the 
construction of our NNLOPS event generator utilize the exact, leading order, 2 —> 2 matrix 
elements. Since the qq —>• Hg matrix element contains no enhanced terms proportional to 
1/Pt, it cannot give rise to corrections to the Sudakov form factor, hence it does not feature 
here. 

The pure 6-quark loop mediated contributions to the pt spectrum are greatly sup¬ 
pressed by a factor m 2 /m? H relative to those owing to the interference of the t- and 6-loop 
mediated amplitudes, hence, our first approximation has been to neglect the former. Our 
second simplification is to use the relatively simpler small quark-mass limit of the scalar 
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loop integrals given in ref. [29] for the 6-loop amplitudes, while using the large-m-i limit for 
the top-loop amplitudes. The small quark-mass limit of ref. [29] is defined by taking m b as 
being small compared to all other kinematic invariants in the process. 

We have validated the following approximate matrix elements numerically in various 
ways. To begin with we have compared them point-wise in phase space to the exact matrix 
elements. In the latter case the evaluation of the integrand in the correction to the Sudakov 
exponent, modulo convolution with the PDFs, with the approximate matrix elements agrees 
with its exact analogue to well within 20% for 20 < py < 150 GeV. For qg —?> Hg and gq —> 
Hq channels the approximate matrix elements are particularly effective, and agreement in 
that case is better than 10% for 10 <C px *C 300 GeV. Since all of our matrix element 
approximations involve the small-m;, limit, they can all be expected to breakdown for 
Pt —> rn b , however, to perform the proposed correction to the Sudakov form factor, which 
involves the cumulant in px starting at its maximum value, this breakdown is unimportant. 
We have also tested the approximate matrix elements against their exact counterparts in 
the evaluation of the full correction to the Sudakov form factor itself, A mf) ( px , Vh )• Due to 
the fact this quantity involves an integral over px up to its maximum attainable value, and 
the convolution of the matrix elements with the PDFs — marginalising the contribution of 
the high px > mt region, in which the gg —>• Hg approximation loses quality — we find 
that results for A mi) {px, Uh ) obtained using the approximate matrix elements are almost 
indistinguishable from those found using the exact matrix elements. 


A. 2.1 Preliminaries 


The exact spin-colour averaged Born amplitude squared, including finite top- and bottom- 
quark-mass effects, M .\and its large-m*, limit, M^^ H , are given by, 


[oo,0] = Gpmjj a 2 s 
H 288\/2 7T 2 


and 


K .[m t ,m b } 
J l 99^H 


T M 


[OG,0] 
gg^H i 


(A.33) 


where, 


T = 


J2 ^2 ( r 9 + ( T <1 - !) / ( T ?)) 

q 9 


arcsm ,/r„ 


/ (t,) = 


Tq < 1 , 
Tq > 1 , 

(A.34) 

with T q = ^Tnjj/mq and m. q being the mass of the quark (q = t, b). In the large-m£-small-m£> 
limit, 


1 l + v/l — 1/ T q 

log - V , - ITT 

6 1-^1-1/t, 


lim T ->• 1 + j r b 1 T t b + O (r b 2 ) , T t b = 4 + 7r 2 - log 2 (4 r b ) . (A.35) 

mt-Kx) 4 

m^—^O 


In eq. (A.35) the leading term (1) arises from pure top-loop mediated gg —>• H amplitudes, 
the second term proportional to owes to the interference of top- with bottom-quark 
loop gg —>• H amplitudes. To evaluate finite mass corrections to the Sudakov form factor 
we neglect terms O (t^ 2 ), arising from pure 6-loop mediated contributions. 

For the pole 6-quark-mass, m b = 4.75 GeV, and mn = 125 GeV we find \J~tb T b l = 
—0.125, while taking instead the MS 6-quark-mass, m b = 3.38 GeV, leads to \Ftb T b l = 
-0.084. 
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A. 2.2 qg —>• Hq and gq —> Hq approximate matrix elements 

Neglecting bottom quark contributions, in the larg e-nit limit the qg 
is given by, 


M 


[oo,0] 

qg^Hq 


= M 


[oo,0] 

gg^H 


SttooCti 


m 


H 


t 2 + s 2 


U 


Hq matrix element 


(A.36) 


Our approximate 2 —> 2 matrix element for qg —> Hq is obtained by taking the large- 
TBi-small-jnf, limit of the full qg Hq matrix element including top- and bottom-quark 
mediated contributions. We express it in terms of a piece, M qg ^ Hq , which respects con¬ 
ventional soft and collinear factorisation, and a further remainder term, which 

does not: 


lim M [ ZT^ a 

m^O 


M 


A 

qg^Hq 


= M 


T 

qg^Hq 


+ M 


'll 

qg^Hq 


where, 


M 


T 

qg^Hq 


= wit 00 ’ 0 ! 


qg^Hq 




(A.37) 


(A.38) 


M 


n 

qg^Hq 


= M 


[oo ,0] 
qg^Hq 


m 


H 


+ m 2 H 


3 r -l 
4 b 


U 


-o“ Ftb + log 


ml 


m 


H 


\u\ + £ 


+ m? H 


, m 

lQ g ~r 


m 


H J 


(A-39) 


with e = m? inserted by-hand, to regulate the potential divergence in the region where our 
small-rrife approximation breaks down. 

We note how, in the approximate matrix element here, the interference of the top- 
and bottom-quark loop contribution vanishes at high |u|. Since the top-loop amplitude 
contribution to qg -A- Hq is manifestly real in the large-rrq limit, the vanishing interference 
at high |w| implies the bottom-quark loop amplitude is pure-imaginary, i.e. the latter is 
dominated by real rescattering at high energy. 

In the limit of large u (backward-scattering of the incoming quark) the large -mt limit, 
severely fails to describe the Higgs boson’s transverse momentum for intermediate pr- On 
the other hand, the quark-mass dependence of the purely top-quark mediated contribution 
can be captured by a relatively simple formula: 


M 


[m t , 0] 
qg^Hq 


M 


[oo,0] 

qg^Hq 


X 


36 mj 

( m H + M ) 4 






The approximate matrix element implemented is given by the rescaling 


M qg^Hq 


M 


[mt , 0 ] 


qg^Hg 

A/ft 00 ’ 0 ! 

JVl qg->Hq 


X M^g^Hq 


(A.41) 
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The analogous gq —y Hq matrix element is readily obtained by substituting u +-y t in 
the qg —y Hq formulae above. 


A.2.3 gg —y Hg approximate matrix element 


Neglecting bottom quark contributions, in the large-?Bt limit the gg —y Hg matrix element 


is given by, 


M 


[o°,0] 
99^ Hg 


M 


[°o,0] 
99^H 


( 87rn s 2 N C \ 

V mi ) 


m H + S 4 + t 4 +U A 
2 stu 


(A.42) 


Omitting terms O (r^ -2 ), in the large-mt-small-m& approximation, we go on and take the 
limit mb Pt ■C m#, keeping only the resulting leading log 2 m 2 term in the top-bottom 
interference piece we obtain 


J^oo M gg^Hg M 99^Hg + M fg^Hg , ( A ' 43 ) 

mj->0 

m b <^p T <^m H 


where, 


M M^Hg - 


M 


[oo,0] 

gg^Hg 


1 + l Ftb 



(A.44) 


M 


n 

gg^-Hg 


M 


[oo,0] 

gg^-H 


87tq! S 2 N c 
Pt 


1 3 

?4 T| 


-i 


(1-z) 




. (A.45) 


For z —> 1 this approximation recovers the results of conventional soft eikonal factorisa¬ 
tion at the level of the pure t-loop mediated contribution, as well as for those terms above 
originating from t-b interference. The vanishing of the factorization-breaking remainder 
pieces as z — y 1, is consistent with Sec. 3.1 of ref. [40] concerning soft factorisation. The 
formulae here extend those of ref. [40], giving a unified description of soft and collinear 
regions. 

Starting from the large-mf-small-mf, approximation and taking instead just the z — y 0 
high energy limit, neglecting terms 0(1/z) gives formulae with a structure strongly re¬ 
sembling of the large-JUf-small-JUf, approximation applied to the qg —y Hq and gq —y Hq 
channels eqs. (A.38, A. 39): 


lim M 

mt—>oo 


[mt,m b ] 

gg^Hg 


0 

z—y 0 


M Tg^Hg + M w^Hg 


(A.46) 


where, 


M 'gg^Hg 


M 


[°o,0] 

gg^Hg 




(A.47) 


M 


in 

gg^Hg 


.,[00,0] rn 2 H 
gg^Hg ^ + 




Htb + log 2 



4 Pt ^ 

- 5 —. -2“ lo g 

+ m A H 


Pt \ 
m 2 H ) 


(A.48) 
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Our final approximate matrix element expression for gg —» Hg is made by combining 
those for the mb <C Pt *C m h and z 0 domains into a single formula respecting both 
limits and interpolating between them: 


a Amt,mb] 

JVl gg^Hg 


M tg^Hg 


M 

M 


[m t , 0 ] 
gg^H 
[oo,0] 
gg^H 


[M^ Hg + Mfg^Hg] 


with 

= i + |w) 


(A.49) 


(A.50) 


M 


n 

gg^Hg 


= M 


[oo,0] 


m 


H 


gg^Hg p 2 + m 


H 


3 r -i 
4 b 


Pt 


m 


J^tb + 


H 


\{\-z + \z 2 ) 

z(l-z) + ^ + ^i 


log 2 


mr 


Pf + £ 


4 p\ 


(A.51) 


1 ( P T 

2 i 2 ( 2 

Pt + fflj \ 


As in eq. (A.39), e = has been inserted by-hand in eq. (A.51) to regulate the spuri¬ 
ous divergence in the region where our small-m;, approximation breaks down. The ratio 
^gg-^H appears in eq. (A.49) as an economical bid to capture some of the finite 
top-quark mass dependence. 
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